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Abstract 

We construct various exact analytical solutions of the SO{‘i) BMN matrix model that correspond 
to rotating fuzzy spheres and rotating fuzzy tori. These are also solutions of Yang Mills theory 
compactified on a sphere times time and they are also translationally invariant solutions of the 
N = 1* field theory with a non-trivial charge density. The solutions we construct have a Zisf 
symmetry, where N is the rank of the matrices. After an appropriate ansatz, we reduce the 
problem to solving a set of polynomial equations in 2N real variables. These equations have a 
discrete set of solutions for each value of the angular momentum. We study the phase structure of 
the solutions for various values of N. Also the continuum limit where Y —)• oo, where the problem 
reduces to finding periodic solutions of a set of coupled differential equations. We also study the 
topology change transition from the sphere to the torus. 
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I. INTRODUCTION 


The BMN matrix model [T] is a massive deformation of the BFSS matrix model [2]. The 
BFSS matrix model describes the discrete lighten quantization of M-theory on flat space. 
The BMN matrix model analogously describes the discrete light cone quantization of M- 
theory on a maximally supersymmetric plane wave. These maximally supersymmetric plane 
wave geometries were constructed by taking a Penrose limit of supersymmetric AdS x S 
spaces in [3]. 

The BFSS matrix model results from the dimensional reduction of A/” = 1 SYM from ten 
dimensions downs to 0 + 1 dimensions and corresponds to the dynamics of DO-branes |1] . It 
has an SO{9) symmetry of the transverse directions and a gauged U{N) symmetry, where 
N is the rank of the matrices. The BMN matrix model splits these transverse directions into 
two sets with different masses, so that the bosonic symmetry reduces to an SO{3) x 50(6) 
subgroup of 50(9). The subset of the theory where only the 50(3) charged scalars are 
excited is the 50(3) BMN matrix model. The full model also results from considering an 
SU{2)l invariant set of conhgurations in Yang-Mills theory on an 5^ x M geometry [5]. 
Thus, any classical solution of the 50(3) BMN matrix model is also a classical solution of 
Yang Mills theory on a sphere times time. One can also show that any solution of the BMN 
matrix model corresponds to a classical solution of the TV = 1* field theory ( see [B] and 
references therein). 

It is known that for generic initial conditions, the BMN matrix model is chaotic [M]. 
This can be understood from the chaotic dynamics of dimensionally reduced Yang-Mills [niT 
[13]. However, special initial conditions can in principle be soluble analytically. It is expected 
quite generally that solutions that minimize the energy with an additional conserved quantity 
turned on and constrained can be stationary. This is usually handled with a Routhian if the 
conjugate variable to the conserved quantity can be separated. Many of these states that 
minimize the energy given some conserved central charge have interpretations in terms of 
BPS states in supersymmetric held theories. 

As it turns out, the BMN matrix model has exact, supersymmetric solutions with zero 
energy [1]. These matrix conhgurations are characterized by all adjoint representations of 
su(2). They have an interpretation as giant gravitons [Tl|. The spectrum of huctuations 
around these solutions is known [15] (see also [16] for an alternative derivation of the spec- 
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trum) and one can argue that there is a large tower of protected states that are available to 
study HT]. Unfortunately the nonlinear structure of the classical solutions that make this 
tower of BPS states is not known. 

It is expected that adding angular momentum to the fuzzy sphere states can induce 
topology changes from a sphere to a torus [IH]. Our purpose in this paper is to investigate 
this topology transition with a special family of matrix solutions at hnite angular momentum. 
The paper is mostly devoted to constructing these solutions. Once the solutions are found, 
the geometry of the corresponding fuzzy membrane is analyzed using the techniques in (TH] . 


II. THE HAMILTONIAN AND THE ANSATZ 


The SO{3) BMN matrix model is a dynamical system with three N x N Hermitian 
matrices or alternatively X, Y, Z. The conjugate momentum matrices are Pi, 2 , 3 , and 

Px,Y,z respectively. The BMN Hamiltonian is given by 

H = + />! + Pi) + ill' j (1) 

The system possesses a U{N) gauge symmetry where X* and Pi both transform in the 
adjoint, X* —)• UX^U~^ and Pi —)■ UPiU~^. The presentation of the Hamiltonian Q is in 
the gauge Hq = 0. The generators of gauge transformations are the matrix of functions on 
phase space given by 
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® = p) 

i=i 

The dynamics need to be supplemented by the Gauss’ law constraint 0 = 0. The system 
also enjoys an SO{3) symmetry of rotations of the matrices X, Y, Z into each other. The 
generator of angular momentum along the Z direction is 


J = Lz = Tr(XPy - YPx) 


( 3 ) 


with similar expressions for the other two S'0(3) generators. Lastly, the equations of motion 
are 


X^' 

Pj 


dH 
dH 




= -X^ - 3iP^^X^X^ 


[[X^X”^],X”*] 
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( 4 ) 

( 5 ) 




The solutions with H = 0 are given by fuzzy spheres. These are solutions of the equations 

[X\ X^'] = (6) 

The solutions to these equations are characterized by direct sums of the adjoint matrices for 
irreducible representations of su(2). For these solutions we have Pi, P 2 , P 3 = 0 and thus are 
classically gauge invariant according to (|^. These solutions carry no angular momentum as 
L = 0 identically. 

The Hamiltonian also satisfies a BPS inequality bound, where H > \ J\ (details can be 
found in [23]). Solutions that saturate the bound will be called extremal or BPS. This 
follows from writing the Hamiltonian as a sum of squares in a slightly different way 

// = IV (^ip| + i(Pi ± (X, + '«23 i1X»,X']))2 

+ ^(P2 =F (A'l + iei23[.V, X=>]))" + i(X3 + fc3i2[X', .V])d ± J (7) 

The cross terms between X 2 , Pi and Xi, P 2 in the squares generate a copy of J that needs 
to be subtracted. The cross terms with Pi and [X^,X^] lead to something that does not 
automatically cancel for generic matrices, but after a bit of reshuffling can be shown to be 
proportional to 

X3([X\P'] + [X2,p2]) (8) 

and we recognize the Gauss’ law constraint starting to arise. After imposing the full Gauss’ 
law constraint, we get Tr(X^[X^, P^]) that does vanish identically. 

The BPS bound is not directly related to supersymmetry. Instead it is derived from the 
conformal group in four dimensions, where one can show based on unitarity arguments that 
E > J from requiring that K = P\ where P, K are the generators of translations and 
special conformal transformations on x M. This bound asserts that the dimension of 
operators is greater than the spin and is usually saturated for free theories, or nearly free 
theories at leading order in perturbation theory. This should be lifted in general theories (see 
[21| for a recent discussion). However, the bound does show up in studying supersymmetric 
BPS states [25l |26| for the full BMN matrix model. This bound descends to any classical 
solution of Yang-Mills theory which is a conformal field theory at the classical level. 

Now we move on to construct solutions to (|^ and ([^ with non-zero J . We dehne the 
matrices 

X* = \(X + iY), X- = \(X-iY) (9) 
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X=X+ + X-, Y = -i{X+-X-) 


(10) 



\ 


( 11 ) 


0 


0 aAr_i exp(ia;Ar_it) 


yoAT exp(za;7vt) 0 


0 


with Oj constants, and X (t) = (X+(t))'^ (this is the transpose complex conjngate). At the 
same time we also take 


Z{t) = diag(2;i, ...,zn) 


( 12 ) 


independent of time and real. We will explain the origin of this ansatz in the following 
section. Notice that all fnzzy sphere gronnd states are solntions of this kind already, with 
Ui = 0, and some of the a* = 0. Gange transformations that commnte with Z allow ns to 
vary the relative phases of the a*. Therefore we can assnme that they are real, and a common 
phase can be translated away by choosing the starting time appropriately. This ansatz is 
different to those that have been stndied before [23] (previons works also look for solntions of 
the BMN and BFSS matrix model where more than three matrices are oscillating [2HI1^). 
Other time dependent solntions can be fonnd in [301 ?' , |3T] , bnt again, these involve more 
matrices being tnrned on. In particnlar, we allow for mnltiple freqnencies to arise in the 
ansatz, rather than jnst one. This does not affect fact that the system is rigidly rotating. 
One can nse a gange transformation to make the freqnencies the same, bnt one pays the 
price that Aq, the connection in the time direction, becomes non-trivial. This is actnally 
very nsefnl for the BPS states, where Aq oc z 125 ]. 

The method of solving the eqnations is then to hrst solve for the conjngate momenta by 
nsing eqnation (|^. Then, we solve the Ganss’ law constraint (|^ to relate the oji to each 
other. These can then be snbstitnted into the angnlar momentnm eqnation ([^, so that we 
can express the Ui in terms of the total angnlar momentnm = J and the a*. One can 
then show that the system of eqnations of motion (|^ rednces consistently to an algebraic 
set of real eqnations for the a*, Zi and that it has as many nnknowns as there are variables. 
One therefore expects to generally find a discrete (possibly empty) set of solntions to the 
eqnations. We will show eventnally that this set of solntions is non-empty for all J. 
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A. Relation to SYM and Af = 1 


The BMN hamiltonian also arises from sphere reductions of four dimensional Yang-Mills 
into a X M and from the Af = 1* theory. The relation to Yang Mills on x M is as 
follows. Consider that the round three sphere is also the group manifold of SU{2), and has 
an SU{2) x SU{2) isometry group by acting with the group on the left and on the right. 
We can use a basis of left invariant one-forms under SU{2), and write the spatial part 
of the connection connection as follows 

A = Aie^ (13) 

where the A^ are now Lie-algebra valued functions on S^. Requiring that the allowed con¬ 
figurations are invariant under left actions of the group, we have that Ai becomes position 
independent and is just a constant hermitian matrix. The Mauer-Cartan equations for the 
one forms e* then give us that 

dA + A f\ A = Aide^ -f [Aj, A (14) 

whereas the electric fields will give 

DtA ~ {DtAfjdt A e* (15) 

Because of the large amount of symmetry preserved, we are led to a consistent truncation 
of the SYM lagrangian. After plugging in, we identify Aj ~ Xj and P* ~ DtAi. The 
Legendre transform of the Lagrangian for SU (2) invariant fields will give rise to the same 
BMN Hamiltonian. 

For the Af = 1* held theories, they are obtained from Af = 4 SYM after an 50(3) 
invariant mass deformation of the super-potential. The Af = 4: SYM Lagrangian can be 
characterized by having three chiral matter super helds in the adjoint of SU{N), and 
a super potential of the form 

W = Tr(0VV^ - 0^0^0^) (16) 

and the full potential has an 50(6) R-symmetry invariance, of which only an SU{3) rotation 
is apparent in terms of AT = 1 super helds. 

TheA/" = 1* deformation adds the super potential mass term 

W = -iMTr(|^(y)A (17) 
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The factor of i is a choice of convention, as the phase of M can be changed by a global R- 
charge rotation of the 0. This preserves an 50(3) global symmetry from the R-charge, but 
the theory can confine in the infrared in some of it’s vacua. For this theory, the classical vacua 
are given by fuzzy spheres, or su( 2 ) representations 120 ] (see also 1211 for the characterization 
of the vacua at strong coupling). The potential for the 0 helds is given by 



( 18 ) 


where the first term is the F-term and the last term comes from the D-terms of YM. The 
theory also has a parity transformation that sends the super held 0* —>■ This is a 
symmetry of = 4 SYM if the theta angle of the held theory vanishes. The scalar potential 


is invariant under this transformation. We can choose to look for conhgurations that are 
classically parity invariant. In that case, we only preserve the real part of (/>*, while the 
imaginary part is removed. This parity transformation commutes with the 50(3) symmetry 
group, but not the full 5f/(3). 

With this constraint on the helds, the potential term arising from the D-terms automati¬ 
cally vanishes. Moreover, the potential for V (0) becomes the BMN potential after factoring 
out the dimensionful constant from the helds 0. A translation invariant classical solution 
of the W = 1* held theory with parity invariance can be understood as a classical solution 
of the 50(3) BMN matrix model. The angular momentum of the BMN matrix model solu¬ 
tions becomes a charge density for a global symmetry of the held theory. This is a charge 
density in a four dimensional A/" = 1 held theory, so it is not a charge density for a central 
charge. The solutions are to be regarded as non-supersymmetric, but in the classical limit 
they are controlled by the same dynamics as the 50(3) BMN matrix model. The quantum 
corrections will be diherent. The structure of solutions will then be described by the (par¬ 
ity invariant) phase diagram for the weakly coupled A/" = 1* held theory at hnite charge 
density. This in turn can be understood as a phase diagram for a (top-down) holographic 
superconductor [ 22 ] . 

III. SYMMETRY CONSIDERATIONS 

A rather natural question to ask is if we can hnd rotationally invariant conhgurations of 
the 50(3) BMN matrix model around the Z axis that rotate uniformly without changing the 
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shape of the conhguration and are a small perturbation of a single fuzzy sphere conhguration. 
After all, this is how the giant torus conhgurations in supergravity are constructed [TH]. It 
turns out that the answer is no. 

The way to see this is the following. Assume that Z is a Hermitian matrix with eigenvalues 
that are not degenerate. A conhguration will be rotationally invariant around the Z-axis if 
a naive rotation of the matrices into each other can be undone with a gauge transformation. 
This is the same way that the method of images works for D-branes on orbifolds when 
considering a discrete subset of the rotation group [31] (this is also the mechanism for 
rotational invariance of monopole solutions in nonabelian gauge theories [331 [33] h 
That is, if for any angle 0 we can hnd a unitary matrix [/(6) such that 

u{e)zu-\e) = z 

U{e)XU-\e) =Xcos(0) - Ysin(0) (19) 

[/(0)YU-^(0) = Ycos(0) + Xsin(0) 

Because Z has non-degenerate eigenvalues, the hrst equation tells us that f/ must be diagonal 
in the same basis that Z is (they commute with each other). We gauge transform to such 
a basis without loss of generality. The angle 0 will be identihed with the time evolution 

itself later on. The uniform rotation motion of the conhguration requires that Z is time 

independent, so that Z = 0. 

Let |1),..., \N) denote the eigenvectors of Z. It is also convenient to rewrite the last two 
equations of ( [I^ in terms of the general matrix X~^ = X + iY and it’s adjoint. The matrix 
X~^ is unconstrained. Then we have that 

u{e)x+u-\e) = exp{ie)x+ ( 20 ) 

Using a general expression for X~^ \m) = |n), and U ~ diag(exp(i6'i)) we hnd that 

transforms under conjugation by U as 

exp(i(0™ - 9n))X+^ (21) 

which can only be equal to exp{i9)Xmn H 9m — 9n = 9 mod (27r) or if X^^ = 0. Picking 9 
arbitrarily small, we can only satisfy 9m — 9n = 9 for some of the components, and because 
9 is very small, we can arrange the kets \n) such that the On are strictly decreasing as we 
increase n. This shows that U{1) invariance requires to be upper triangular with zeros 


on the diagonal. Now, comparing with the fuzzy sphere solutions, X~^ actually must border 
the diagonal. That is, we hnd that the solutions must be of the form 




^0 |ai| exp(i0i) 0 ... ^ 

0 0 |a2|exp(z02) 

0 ... 0 |aAr_i|exp(i07v-i) 

\ 0 0 ... 0 / 


( 22 ) 


where the Oj are some (as of yet unspecihed numbers, that depend on time), and both X~^X~ 
and X~X~^ are diagonal and rotationally invariant themselves. Thus if a conhguration is 
rotating uniformly we hnd that the only possible solution has the |aj| being constant, and 
all the dynamics will be in the phases 0i(f). Inserting these expressions in the Gauss’ law 
constraints, we hnd that 

loilVi - ki+ilVi+i = 0 (23) 


for all i. This leads us to show that 4>n-i = 0 and from there, (pi = 0 for all i. The 
solutions are thus static and carry no angular momentum. Since we want to turn on angular 
momentum, we need to relax the rotational invariance around the Z-axis of the conhguration. 

The obvious idea is to look for the maximal discrete subgroup of rotations that can 
actually be preserved if we can not have a full SO{2) symmetry. Let us assume that we turn 
on an X^n which is not one of the above. The matrix U that implements the constraint 


(20) for X^ dehned in (22) is diagonal and up to a global phase it is equal to 


U = diag[exp(—i6*), exp(—2i6*), exp(—3i6*)... exp(—A^i6*)] (24) 

We then have that nO — mO = 9 mod (27r), or equivalently, that (m + 1 — n)9 = 0 
mod (27r). This tells us that we can preserve a Z^+i-n subgroup of the S'0(2) rotations if 
we turn on this particular X^n- If we want to maximize this group, we hnd that we must 
take m = N and u = 1, where we get an Z^r subgroup of the rotation group to be invariant. 
None of the other Xmn that are not already turned on are neutral under this subgroup. We 
can self-consistently set them to zero by symmetry arguments; if an initial solution for X, P 
respects the symmetry, the dynamics being Z^v invariant will guarantee that no symmetry 
breaking can occur afterwards. We are therefore led to a general ansatz for X+ which can 
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still have arbitrary time dependence 




0 

0 

0 

Oat I exp(i0Ar) 


I Oil exp(i0i) 
0 


0 

|a2|exp(i02) ••• 


0 

0 




0 


0 |aAr_i| exp(i0Ar_i) 
0 


(25) 


Now we want to insist that the only motion that the system does is rigid rotation, so that 
the \ai\ are necessarily constant, but the can vary in time. Substituting this in Gauss’ law 
shows that all the motions in the angles are related to each other by — laj+ip^j+i = 0 

cyclically. Furthermore, one can show that ~ J which is a conserved 

quantity. Putting these two pieces of information together shows that the phases have 
constant time derivatives which we call cuj. This leads us to the general form of the ansatz 
described in equations ( pT| and (12). 

Notice also that the ansatz we have made has an additional remnant discrete Ztv symme¬ 
try on the variables a^, Zi, where we send |afc|exp(i0fc) —)■ |aA;+i| exp(i0fc+i) and —)■ Zk+i 
cyclically. That is, any solution of the equations can be permuted to a new solution using 
this symmetry. This is a subset of the gauge freedom of the original system, where we 
permute the eigenvalues of a hermitian matrix. In DO brane dynamics this permutation 
symmetry is associated to the permutation statistics of D-branes [2]. This symmetry can 
also be understood if we orbifold the matrix problem by the original symmetry we iden- 
tihed via the rules of [M]. In the orbifold theory by an abelian symmetry, it is expected 
that one has a dual quantum symmetry that permutes the nodes of the corresponding 
quiver theory EH. Gauging this dual quantum symmetry restores the original theory in a 
straightforward way [3H]. There is a second symmetry where we reverse the order of the 
|ai|, Zi and also change the sign of the Zi. This acts essentially as reflection on the Z-axis, so 
that when we combine it with time reversal, we can still spin in the same direction. These 
two symmetries form a dihedral group with 2N elements and it will be useful for analyzing 
the set of solutions of the ansatz. 

Notice that if a discrete subgroup of the quantum symmetry is left unbroken, this implies 
that there is an enhanced unbroken gauge group so long as it is not just the spatial reflection 
symmetry. For example, if the Z^r is unbroken (we hnd some solutions of this type), the 


10 





unbroken gauge group turns out to be The solution can be interpreted as N DO 

branes separated from each other in a Zjv symmetric pattern of rotations around the origin, 
just like one would expect from the method of images in an orbifold. 


IV. THE SOLUTIONS OF THE ANSATZ ARE A SET OF CRITICAL POINTS 
OF AN ENERGY FUNCTION 

The equations of motion that follow from the Hamiltonian Q come in two different sets. 
First, we have that 

= Pj (26) 

We can solve these immediately given the ansatz for the X. We hnd that Pz = 0 and that 

Px = X = iVLX+ - iX-VL (27) 

Py = Y = DA+ + X-VL (28) 

where the matrix O of angular velocities is given by 

D = diag((:ni, c<;2, • • •, ojn) (29) 

This way we hnd that 

PxX - XPx = (fDX+ - fA-D)(X+ + X-) - (A+ + X-){iVtX^ - iX-^) (30) 

PyY - YPy = (DA+ + A-D)(-i(X+ - V)) - (-f(A+ - X-))(DX+ + X-D) (31) 

Adding these two, we hnd that the terms with two X+ or two X~ cancel each other, so that 
the Gauss’ law constraint reads 

Q = 2iVLX^X- - AiX-VtX^ + 2iX^X-Vt (32) 

A straightforward computation shows that this is a diagonal matrix and that it is equal to 

Q = 4zdiag(|aj|V - |ai+i|V+i) (33) 

To satisfy the Gauss law constraint, we need that Q = 0 identically. Hence we hnd that 
lajpcuj is independent of i. 

Similarly, we hnd that the angular momentum is written as 

J = Tr((X+ + X-)(DX+ + VG) - (-iX+ + iX-){iVtX^ - iX-Vi)) (34) 
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and similarly, the terms with two copies of X~^ or X cancel even before taking the trace. 
We hnd that 


J = 2Tr(X+X-fi + = 2Tr(diag(a;i|ai|2 + (35) 

Now, notice that becanse of Ganss’ law constraint, all the u:i\ai\^ are eqnal to each other. 
Thns, even before taking the trace, the matrix version of J is proportional to the identity. 
We can interpret this as having nniform density of angnlar momentnm per nnit DO-brane 
of the corresponding fnzzy membrane. We nse this to hnd that 



J = 4:NuJi\ai\'^ 

(36) 

so that we can substitute 



J 

~ 4N\ai\^ 

(37) 

Knowing the Oj , the Zi 

and J, we can evaluate the energy by substituting the above results 

in the Hamiltonian ([^. 

The result for the kinetic energy is 



l«il) = + Py) = + tPy)(PY - tPr) 

(38) 


= ^Tr(GX+X-G) 

(39) 


N 1 ^ 72 

2=1 2=1 

(40) 


What is important to realize is that this looks very similar to an angnlar momentnm cen- 
trifngal potential, where each particle (associated to the radial variable a*) has the same 
angnlar momentnm J/N and the same mass (in this case the mass wonld be interpreted 
as 4). To specify the fnll problem, we need to evalnate the potential energy as well. A 
straightforward, thongh rather tedions procednre shows that 

^ 1 

G(|aj|, Zi) = -[zj + 2|ai_ip — 2|aj|^]^ + 2(1 + Zj+i — ZiY\ai\^ (41) 

i=l 

where the i are dehned modnlo N. We denote the fnll energy of a conhgnration (J, |aj|, Zi) 
by 

T(T, IcijI, Zj) Ekin(yJj II) T G(IOjI, Zj) (42) 

When we consider the energy fnnction S{J, |aj|, Zi) and recall that = 0 in our ansatz, 
it is straightforward to notice that the equations of motion of the Z variables are exactly the 
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equations that extremize as a function of the Zi keeping the other variables hxed. Namely, 
that dz£ = 0. This suggests that to look for solutions of the original problem we set out to 
do with our ansatz, it is enough to consider the extrema of the £ energy function at fixed J. 
Indeed, any solution of the ansatz that solves the equations of motion of the Hamiltonian 
Q are going to be extrema of the energy function £ and vice versa, any extremum of the 
energy function can be shown to give a solution of the equations of motion derived from the 
50(3) BMN Hamiltonian. 

The essence of the proof is that when we take the derivatives of the kinetic energy with 
respect to the |aj| we hnd that 


d\ai\Ekin 




Am\ai 


= —AujlaA 


(43) 


and these can be assembled into Px = X, Py = Y. Then we need to compare this expression 
to the derivatives of the potential 

( 44 ) 


which can be assembled into the right hand side of the Hamilton’s equations for ([^. The 
two sets of equations can be shown to be the same set when we remove the time dependent 
phases exp(za;if). 

The potential is a sum of ‘nearest neighbor’ terms and it is a sum of squares. As a 
function of the Zi it is quadratic. Therefore if we fix the \ai\, we can solve for the Zi via a 
linear set of equations, and these equations are independent of J. The quadratic form that 
appears in front of the Zi is of the form 6 ij + Uij where H is a non-negative matrix. Therefore 
it is always invertible. 

The kinetic energy is also a sum of squares. Therefore the energy function is bounded 
from below. Notice that when J 7 ^ 0, and as we take \ai\ —)• 0 the energy diverges at |aj| = 0. 
Also, when we take the |ai| —)■ 00 the potential grows quadratically in the |ai|, and if we 
solve for the Zi in this limit, they are bounded. Therefore we expect that generically the 
potential grows quartically at inhnity. The only time when this does not happen is when 
we take a double scaling limit where |ajp = A —>■ cx) independent of i. In this case one has 
that the —)■ 0 and the potential grows only like A^. 

Considering the conhguration space as the set of open intervals |ai| G (0, cxd) (since after 
all we can solve for the Zi given the Oj), we have that the domain of interest is an open 
ball (it is diffeomorphic to ( 0 , 1 )^), and the potential function diverges on all the boundary. 
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while it is finite (and indeed analytic) in the interior. This shows that the energy function 
has at least one minimum. Moreover, if we compactify the boundary of the configuration 
space by adding one point at infinity, we get a configuration space which is compact and has 
the topology of a sphere. The point at infinity realizes the maximum of the energy function 
continuously 

Notice that 


dz 

I) 1-2 -^2 (^2) 7 ^2) | ^2’ ^0 L 


(45) 


since the Zi solve the = 0 equations. Any critical point of the original 

£{J,ai,Zi) will give rise to a critical point of £j{\ai\) and vice versa. In our original ansatz 
where we have ?)N variables, given by ai^Ui, Zi, we have managed to reduce the problem to 
a set of algebraic equations in N variables, the themselves. 

For the purpose of analysis of configurations, the function £j{\ai\) where we fix the J 
and have already solved for the Zi will be thought of as a Morse function on this topological 
sphere (the reader unfamiliar with Morse theory should look at [351 SO] , and the lecture notes 
by Hutchins are very approachable BB)- The main reason for using Morse theory is that 
the main ingredient that it utilizes and analyzes is the set of critical points of a function, 
and from this set the topology of the manifold can be reconstructed. For us, the topology 
is already known (it is a sphere), but the set of critical points is not. The set of critical 
points and gradient flows between them gives a model for the cohomology of the manifold, 
which in our case is known. Given a set of critical points we can ask if the set is consistent 
with the topology of the manifold. If it is not, we are missing critical points. Also, as J 
changes, the set of critical points can change dimension and how these changes can happen 
is understood in general. In our case, because of the extra symmetries of the potential, the 
critical points will exhibit also a representation of the symmetry group. 


^ We can always map the energy function from [0,oo) to [0,1] by using tanh(£) for example. 
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V. THE CASE OF 2 X 2 MATRICES 


The case of iV = 2 is the simplest we can analyze given the structure of our ansatz that 
is not completely trivial. The matrices take the form 


X+ = 


z = 



0 |ai| exp(ici;it) 

|a2| exp(ia;2t) 0 

When we compute the equations that the Zi must satisfy for a*, J hxed, we hnd that 

0 = Zi |ni|^(—6 + 4:Zi — 42 : 2 ) + |n2p(—6 + Azi — 42 : 2 ) 

0 = ^2 + |c['iP(~6 + 4^2 ~ 4^1) + |c['2p(~6 + 4^2 ~ 4^1) 


(46) 


(47) 

(48) 


Summing the two, we hnd that Zi + Z 2 = 0, and then we can substitute this result back, to 
hnd that 


2 : 1,2 = ± 


^ 6(|aip 

) 

02 

P) 

1 + 8 a] 

.P + 8|a2P 


(49) 


At this point it becomes clear that the expressions simplify if we consider the two variables 
P = |aip + |a 2 p and Q = |aip — |a 2 p. This is because P is even with respect to the 
Z 2 symmetry |ai| <H- |a 2 |, and Q is odd. Thus, the symmetry algebra acts simply on the 
variables P and Q themselves. We hnd then that 

J2 32(P - i)g2 


T ( T, (IOj I), I Ui I) 


+ 


+ 


+ 2P 


(50) 


16P - 16Q 16(P + Q) ' 8P+ 1 

When we compute the equations that P, Q must satisfy, obtained by considering dp^qE = 0, 
it is convenient to eliminate the J dependence of one algebraic combination of these. After 
a bit of work, this is accomplished by considering 

+ (P- + Q^)a,E ^ (51) 

Notice that this factorizes, so there are three branches. One where Q = 0 identically, one 
where P = 1/4 identically and another one where 

-32P3 + 4P2 + P 


O = 

^ 16(2P + 1) 

which is positive only if P < 1/4. 

Let US analyze the hrst one. We can substitute Q = 0 in P, to hnd that 

P 

E = — + 2P 
8P 


(52) 


(53) 
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And the minimum occurs for P ^ J/A (here we have taken J > 0, and obviously P is 
positive since it is a sum of squares). We can then evaluate that for this solution 

E = J (54) 


in the limit J —)■ 0, this solution reduces to the trivial solution where all matrices are 
identically zero. For the second solution, we take P = 1/4, and we hnd similarly that 

+ (1 - i6g2)2 


E = 


2 - 32g2 


(55) 


which is of similar form if we use the variable x = 1 — 16g^ (this is, of the form Ax + Bx ^). 
Again, the minimum occurs when 

Q = (56) 

and we also hnd that E = J identically. This is the solution where we choose to take g > 0. 
There is a similar solution with g < 0 that is a Z 2 rehection of this solution. When J —)■ 0, 
this is the standard fuzzy sphere of 2 x 2 matrices. When J —)■ 1, this matches our other 
solution with g = 0. The three solutions meet at J = 1. Beyond J = 1 this solution does 
not exist anymore. 

In the third branch, we have that 


E 


2J2(2P + 1) 

(4P + 1)(16P-1) 


- 4P^ + 9P - 


9P 

2P + 1 


(57) 


and it is easy to solve for J as a function of P (essentially solving dpE = 0), giving us 

(1-P)P(64P2 + 12P-1)^ 

= (2P+l)^(8P+l) - 


(58) 


So that we end up with a parametric solution J{P), rather than the other way around. This 
can be inverted numerically. This only makes sense if .P > 0, so that necessarily P < 1, 
and remember also that 1/4 > P > 0 from the reality of Q. But moreover, |g| < P, and 
this restricts P to be bigger than 1/16. Notice that when J —)■ 0, there are various values 
of P that can arise as roots. The only one that is new corresponds to P = 1/16. This is a 
fuzzy sphere at half radius. This can be easily understood if we make an ansatz of spherical 
symmetry for a saddle point. In this case the matrices X, Y, Z are all proportional to 
the corresponding Pauli matrices with proportionality constant r. Because the energy for a 
static conhguration is quartic in the matrices, and we have that P = 0 at r = 0 and P = 0 
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at r = 1, and always E' > 0, then the energy mnst be proportional to r^(l — r)^. This has 
a maximnm at r = 1/2, which is the sphere at half radius. This solution also matches the 
solution at Q = 0 that we already had when we set P = 1/4. There is similarly a reflected 


solution with Q < 0, where we take the other square root branch cut of equation (52). The 
new solution of the fuzzy sphere at half radius migrates to higher angular momentum as we 
increase P from 1/16 to 1/4 and is a saddle point. Therefore it has Morse index one. This 
solution and the one that is reflected by taking Q —)■ —Q can cancel the two minima from 
the BPS solution as they meet at J = 1 with the trivial solution that has Z = 0 throughout. 

The full set of solutions is depicted in hgure [Tj There we can see that for low angular 
momenta there are hve critical points. Three are minima and two are saddles with Morse 
index one. The saddles and two of the minima are reflected into each other by the Z 2 
symmetry, and one minimum if at the hxed point. The two saddles and two of the minima 
annihilate each other when J = 1. Because the minima that annihilate with the saddles are 
BPS, the saddles need to approach the extremal limit and thus must touch the hxed point 
set, because one can not descend further from the saddle to the hxed point otherwise. 



FIG. 1. Parametric plot of the solutions for oi, 02 derived either from equation (52), from P = 1/4 
or from the solntion Q = 0 with oi = 02 - Superposed we find level sets of the energy function at 
J = 0.3, with the energy levels shown, and we see that the curves pass through the critical points 
of the energy function. 
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VI. THE CASE OF 3 X 3 MATRICES 


The equations that the ai and Zi need to satisfy can be directly derived from (43) and 


(41). We have not been able to solve them algebraically in general, although there is one 


trivial solution with Zi = 0 and all a* equal to each other. This solution exists for any value 
of the angular momentum (and actually all values of N as well). Our strategy for solving 
the problem is to start with known solutions at J = 0 and perturb them numerically slowly 
by varying J until a new solution near the old one is found. Because the solutions are saddle 
points of the energy function, which is considered as a Morse function, small perturbations 
of the function preserve the saddles in general. These are saddles for any J until a subset of 
the saddles collides. A saddle with index m can be annihilated by a saddle of index m + 1 
or m — 1 0 

There are obvious solutions at J = 0. These are the fuzzy sphere vacua. For any such 
fuzzy sphere with k x k matrices, we can also find an unstable fuzzy sphere at half size, 
which is analogous to the one we found for 2x2 matrices. We can then combine these 
solutions together into new solutions. The reason for this is that in any fuzzy sphere for 
k X k matrices we have = 0 and the matrix is upper triangular. This is true regardless 
of if the fuzzy sphere is stable or unstable. We can then mix and match these solutions and 
put them in some order. For 3x3 matrices this does not matter as there are not too many 
ways of partitioning 3 into integers, but the strategy works in general for other values of 
N. The value J = 0 is technically a singularity of the family of morse functions, because 
the fuzzy sphere vacua end up with some a* = 0, and we argued before that these points 
need to be identified with each other in the one point compactification of the open intervals 
|aj| G (0,oo) to get a sphere topology. Numerically, we can start with these solutions and 
perturb entries that start at zero by a small amount, while at the same time turning on 
a small amount of angular momentum. All solutions persist under this procedure, so we 
can safely describe them by taking the limit J —>■ 0. We also need to count them with 
multiplicity, because we have our dihedral group of (quantum) symmetries that let us find 
new solutions of the Oj by permuting them clockwise, or by reflection. This is described in 
table I If we take the solutions as above and we compute the Morse polynomial with just 


^ Recall that the index of a critical point is number of negative eigenvalues of the Hessian at that point. 
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Splitting 

Stability 

Morse Index 

Multiplicity 

Unbroken Quantum Symmetry 

1+1+1 

S, S, S 

0 

1 

Ds 

2+1 

S,S 

0 

3 

1j2 

2+1 

U, S 

1 

3 

Z 2 

3 

S 

0 

3 

Z 2 

3 

U 

2 

3 

Z 2 


TABLE I. Table of sphere solutions. The splitting indicates the size of matrices of the fuzzy spheres, 
S, U indicates if they are stable BMN vacua, or unstable spheres at half size, the Morse index is 
the number of negative modes of the Hessian (after perturbing by a small j, and the multiplicity 
is the number of copies of the solution that are obtained from using the group symmetry actions 
on a given solution. 

these solutions, we find that 

Mtriaiit) = l + 3 + 3t + 3 + 3t^ + t^ = 7+ 3t + 3t^+ t^ (59) 

where the last entry (the one for t^) corresponds to the maximum of the energy function at 
inhnity. The Poincare series of the three sphere is 

P{t) = l + t^ (60) 

The Morse inequalities require that Mt^iaiif) — P(t) = (1 + t)Q(t) where Q should be a 
polynomial with positive integer coefficients. In particular we should have that M(—1) = 
P(—1). This is not the case. This indicates that we are missing critical points of the energy 
function. Since all the solutions above have an enhanced U{1) rotation symmetry, it makes 
sense to look for solutions with such a rotation symmetry for more saddles. We do this 
graphically in hgure In the hgure the Zi have already been solved for, but the |aj| are 
variables. 

We see from the hgure that there are additional saddles that are not rehection symmetric 
with respect to the diagonal, nor are they on the edges of the graph. These new saddles 
have index 1, and the orbit under the symmetry group produces 6 saddles in total, and no 
unbroken symmetry. These would add an additional 6t to Mtriaiit). With this additional 
set of solutions we hnd that 

M(t) = l + f3 + (l + t)(6 + 3t) (61) 
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1.0 



0.0 


FIG. 2. Energy function at J = 0, loal = 0. The axis indicate |ai|, |a 2 |, and the colored contours 
have their energy values indicated. The solid dots indicate approximate positions for the various 
saddles. The solid dots in red indicate new saddles that do not arise from collections of round 
fuzzy spheres. 


and now we satisfy the Morse inequalities. This means that it is consistent if these are all 
the saddle points and we are not missing any. 

We have not found any other saddle point numerically at J = 0. If they exist they should 
have all three |aj| 7 ^ 0. It is possible to use this information to show that no such saddle can 
exist, and not just for N = 3 but for all N. 

The energy function is given by 

J2 ^ 1 N I 

£{J, \ai\, Zi) = 1 ^ |2 + 2 2|aj_ip — 2|ajp)^ + 2(1 + ZiJ^i — ZiY\ai\^ (62) 


The equations of motion for the ai yield 

J2 

0 = —^*+1 ~ (^*+1 ~ + 2(—|ai_ip + 2|ajp — |aj+ip)) (63) 


Next we suppose that |aj| 7 ^ 0 for all i at J 
divide the rest by 4|ai|. We are left with 


0. The hrst term of (63) vanishes and we may 


0 — (1 + Zi^i — ZiY + (^j+1 ~ Zi) + 2(—|aj_i|^ + 2|aj|^ — |aj_|_i|^) (64) 
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Summing over all i we have 


N 

0 = iV + — ZiY (65) 

i=l 

Since the summation is non-negative and N is positive, we have reached a contradiction. 
Thus 10,1 = 0 for at least one z at J = 0. 

The equation of motion for the z, yields. 

0 = Zi- 4|aip Q -h Zi+i - Zi^ + 4|ai_ip Q + (66) 

Summing these equations yields a traceless condition, Zi = 0. This had to be true because 
Pz = 0 and the trace of the matrix model is just a harmonic oscillator. 

This feature has the following implication. When we perturb away from J = 0, the |ajp 
will be modified very slightly (at order J^), but for the one that begins at zero it is 
different. In the kinetic term we will have a singular term, so that the energy function will 
look like 

£(J) ~ T(0) + a\ae\‘^ + g^2|Q^|2 (^7) 

and clearly the \a(\ that minimize this are such that \a(Y — J/VSN'^a. When we plug this 
into the energy function we hnd that there is always a linear term in J. This means that 
the fuzzy conhguration built with out ansatz can never be considered as a rigid body; for 
rigid bodies the energy goes like J^, where the coefficient of proportionality depends on the 
moment of inertia. The fuzzy configurations rotate by turning on wave-like excitations on 
the sphere. This is exactly as expected for a membrane. 

A presentation of the solutions found at J = 0 for 3x3 matrices, followed as we change 
J is shown in figure As shown in the figure, there are many phase transitions. The one 
marked A corresponds to two saddles of index one and one saddle of index 2 merging into 
a single saddle of index 1. This is repeated in three different locations due to the symmetry 
operations. The unbroken symmetry of the incoming saddles of index one changes from 
all symmetry broken to a saddle with a Z 2 symmetry unbroken. It also shows that the 
saddle of index one ends smoothly and without corners. Actually this is what is expected 
in general: saddles should end smoothly and without corners as we pass through critical 
values of J. The phase transition marked as C is inconsistent as shown: a single saddle 
of index zero and a saddle of index one can not annihilate into a saddle of index zero 
(this would violate the Morse inequalities). This means that there is a solution missing. 
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FIG. 3. Phase diagram of solutions as a function of J on the abscissa, and the |oi| are plotted on 
the ordinate. We follow the solutions slightly perturbed from J = 0. Shown are the values of the 
|aj|. Different solutions for |ai|, |a 2 |, jasl are colored differently, but all values are shown in the 
same color for the same solution. The numbers attached to saddles are the Morse index. Phase 
transitions where solutions merge or end are shown as A, B, C, D. The ones that are not marked 
below the maximal symmetry solution are C, D also. 


Similarly, in phase transition D, the line actually ends. The vertical line down is an artifact 
of joining the numerical solutions, and instabilities of the numerical method when the Hessian 
is degenerate. This solution ending also indicates a missing solution, because of inconsistency 
with the Morse inequalities. Finally, at the transition B we see that a set of three saddles 
of index zero and index one merge with the saddle with all symmetries unbroken. Although 
this is in principle allowed, it is not smoothed out like in transition A, or as one would 
expect from transition C once a solution is found. We can conjecture that the solutions 
that correspond to the green and purple lines in the hgure actually continue to the other 
side. Because of the way solutions are merging and the dihedral symmetry being restored, 
transition H is a multi critical point. It is natural to believe that the solutions that are 
continued end on the other inconsistent transitions. They should have also an unbroken Z 2 
symmetry each and come in three copies. Because of the unbroken Z 2 symmetry, two of the 
\ai\ should be equal to each other. This makes it possible to guess the solutions by looking at 
energy contours in a two dimensional plot where we set jasl = |a 2 | a bit after the transition 
at J ~ 2.1. We should be looking for a BD line of index one (it’s the only way we can 
cancel a solution of index zero). The green line actually saturates the BPS inequality, and 
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FIG. 4. Energy contours in |ai|, |a 2 | = josl at J = 2.1. The saddles are hard to see. The energy 
contours have energies starting at E = 2.1 and spaced at 6 E = 2 x 10“^ showing ten contours. 



FIG. 5. Full Phase diagram of solutions as a function of J on the abscissa, and the |aj| of a single 
solution are plotted on the ordinate. 

this should persist in the analytic continuation to the right. This indicates that there should 
be a saddle of index 0 joining BC. This process is depicted in hgure Once we include 
the new saddles after this point, we can complete the phase diagram in hgure What we 
see from the hgure is that hrst, there is a maximal angular momentum after which there 
is only one solution, and this solution preserves the maximal dihedral symmetry. Second, 
the trajectory of the maximal sphere goes through a phase transition at a hnite J before 
reaching this maximal angular momentum. We have checked that after the hrst transition 
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the solution stops being extremal (BPS), even though it is a local minimum. This means 
that the corresponding phase can be considered metastable from energetic considerations. 

At the maximal angular momentum, it becomes classically unstable. The jump from the 
metastable to the stable conhguration is a hrst order transition, but at the place where it 
changes from being BPS to being non-BPS, we have coexistence of a second order phase 
transitions and a hrst order phase transition with the same energy. This degeneracy is 
expected to be lifted by quantum corrections. At small h we expect that the phase diagram 
has the same structure we have shown, because it is controlled by topological aspects of 
a Morse function. This is similar to the transition found in 021 , where studying spherical 
D-brane conhgurations, beyond the maximal giant graviton there is a non-BPS metastable 
D-brane solution that continues for a while with larger R-charge (this branch was hrst 
found in |13], but the stability analysis was not done), and the family of solutions ends in 
a transition that should take us to another branch. This transition should be where the 
metastability is lost. 


VII. OTHER EXAMPLES 


Our next task is to understand how to go beyond N = 3. Armed with the information 
that at least one of the |aj| needs to be zero at J = 0, suggests that we can begin by looking 
at saddles for smaller N and ht them into saddles for N by either bordering by zeros, or by 
taking direct sums of solutions for smaller N such that the rank adds up to N. We choose 
these saddles to have \ak\ =0 for the entry on the hrst row at the bottom corner. These will 
all ht the ansatz where X~^ is upper triangular, bordering the diagonal, and with a zero in 
the bottom corner. This actually always produces correct solutions at J = 0. From equation 


(66), we hnd that for the values of |aj| such that \ai\ = 0, Zi only depends on the previous 


ones, and therefore decouples from Zj+i, and similarly Zi+i will be related to Zi+ 2 , but not 
the other ones. Secondly, the equations for |aj| are trivially satished in this case, because 


of the I a* I appearing in front of it in equation (63). The rest of the equations are satished 


if they were satished for smaller values of N. This gives us an ample trove of solutions 
with which we can explore the phase diagram. This is explicitly shown in hgure We see 
qualitatively the same structure. There is a maximal angular momentum beyond which only 
the trivial solution exists, and it ends in a hrst order phase transition with a metastable 
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FIG. 6. On the left, partial phase diagram of 4 x 4 solutions as a function of J on the abscissa, and 
the I Oil of a single solution plotted on the ordinate in the same color. Notice the similarity of the 
transitions to the ones for 3 x 3 matrices in figure On the right, we evaluate also the Zi values 
for the given solution of the |aj|. The figure is not symmetric under Zi —>■ —Zi which indicates that 
many of the phases are not ‘parity invariant’ with respect to the dihedral discrete symmetry group. 

phase. Because there are many more solutions at J = 0 at each N (the number is bounded 
below at least by twice the partitions of iV), the full phase diagram is more complicated as 
we increase N, and it is clear from the Morse analysis that we are missing quite a number 
of saddles. Although two new solutions always exists for any N^ the maximal fuzzy sphere 
and the sphere at half size, one can expect in general that there are quite a number of 
intrinsically new (indecomposable) solutions that appear at any N. 

For illustration purposes, we also show some of the solutions for 10 x 10 matrices in figure 



FIG. 7. We evaluate Zi{J) for three different solutions of 10 x 10 matrices, the maximal fuzzy 
sphere, a near maximal fuzzy sphere bordered by zero, and a solution for a fuzzy sphere of spin 
7/2 0 1/2. 
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This shows the procedure for building up solutions by bordering by zero or adding up 
previous solutions in more detail. It is clear that the full phase diagram is quite complicated, 
some of the solutions ending in hrst order phase transitions, and some others ending in what 
appear to be multi-critical points. Here we basically show that the problem is amenable 
to computer calculations. We should point out that for larger N we do not have a good 
strategy to search for the intrinsically new solutions yet. The graphical method that worked 
for = 3 is unsuited for higher dimensions, as we can not visualize the data. 


VIII. LARGE N 


An interesting way to proceed to large N is the following. Since the energy function in 


equation (41) is of nearest neighbor type, we can think of it as a discretized version of an 
energy which is an integral of an energy density. Indeed, this is what the interpretation of 
the solutions as of matrix mechanics as discretized membranes indicates we should be doing 
[milS]. With that in mind, we want to replace Zi —>■ z(9) and the same for |ajp —)■ \d\'^{9), 
where 6^ is a periodic coordinate with period one, rather than a discrete set with N elements. 
We basically take i zz NO. Then expressions of nearest neighbor differences get replaced by 
derivatives ZiJ^i — Zi ^ N~^d 0 Z. We also need to remember that the maximal fuzzy spheres 
are of size N when we have N DO-branes, and we want to rescale this out of the energy. 
Therefore, we write Zi ~ Nz{9), and |ajp ~ N‘^\a\‘^{9) and = N J dO. When we do this, 
we hnd that 


Vpot — 


dO 


- (z - 2(9e|ap)^ -h 2 (1 dez) 


2 I |2 

a 


( 68 ) 


and similarly, the terms that contain angular momentum, when we rescale J —)■ give 

us 

Ekin = N^ j (69) 

The point is that in this rescaling, we get a common factor of in front of the energy that 
can be dropped if we want to, to obtain a classical membrane energy which is independent 
of N. The value of N can be changed effectively by changing the periodicity of 9 without 
changing the energy density further. We use the relation Ntot = N J dO to convert the 
changing period of 6 into a different value of N, and we can always change scales by taking 
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9 —)■ a9, z ^ az and |ap —)■ a‘^\a\^. The energy function also has a Z 2 symmetry where we 
change 2 : —)■ — 2 : and 9 —)■ —9. 

The energy function is now a local integral of the functions 2 :, a and their derivatives. 
The condition to be a critical point of the effective energy is a pair of differential equations, 
one for 2 and one for a that come from the variational principle. A complete set of initial 
conditions requires specifying 2 ( 0 ), a^(0), dgz{0), deo^i^O). These can be evolved in 9, but 
the trajectories need to be periodic with a hxed prescribed period (which we are choosing 
to be set equal to one). This effectively quantizes the set of possible solutions so that they 
are discrete. 

Notice that the system is translation invariant in 9, thus there is a trivial one parameter 
family of solutions that is obtained by translation. This becomes a U{1) symmetry that was 
only realized as a Zat quantum symmetry at hnite N, and which is spontaneously broken 
on most of the solutions. Because \a^\ is bounded for interesting solutions, we can always 
start out from a place where (9e|ap = 0 . 

The fuzzy spheres at zero energy appear at J = 0, and are at the zeros of Vpot- These 
occur at 


dez = -l (70) 

de\a\^ = (71) 

so that they satisfy a first order set of equations (typical of BPS states), rather than the 
usual second order equations. At hnite N, the Z matrix for fuzzy spheres is the matrix of 
spins in an N dimensional representation of 5 u( 2 ). Consecutive matrix elements differ by 1 
and so the equation for z{9) is capturing this effect. If we start from 2 = 0 at 6 ^ = 0 (so that 
z{9) = —9), and |ap(0) = it is straightforward to integrate these and hnd that 

de\a\^ = -^0 (72) 

so that 

wm = aL. - 0V4 (73) 

we hit a^(9o) = 0 at some hnite value of 6 * 0 . This is a singular point in the diherential 
equation system because the ehective derivative term squared for 2 vanishes at this point. 
We are allowed to have large jumps in 2 at this point because there is no energy cost to it. 
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We recover this way that the set of possible ground states is a collection of fuzzy spheres. 
Quantization then requires that each of these have an integer amount of DO brane charge. 
This is reminiscent of the LLM droplet picture |1B], where quantization of the area arises 
from a Dirac quantization condition. Notice that the solution can also be written as 

wm = - z\e) (74) 


which gives |ap + |; 2 ^| = const, as one expects from a sphere written in cylindrical coordinates 
(the additional angle is associated to the Zjv —t rotational invariance of the set of solutions 
we are considering in the large N limit). 

At this point, going beyond the BPS solutions, we want to change perspective and think 
of the variable 0 as a time coordinate, and the effective energy function we had before as 
a Lagrangian whose variational principle gives some non-trivial equations of motion. Using 
this change of point of view, we see that the Lagrangian has a (repulsive towards infinity 
and the origin of |ap when j ^ 0) potential of the form 




eff - 




(75) 


there is a non-trivial curved metric 


ds^ = A{d\a\^Y -I- A\a\^dz^ 


(76) 


with translation symmetry in and a non-trivial magnetic potential associated to the one 
form 

A =—2zd\a\‘^ + A\a\^dz (77) 

which in these coordinates produces a constant magnetic field. Notice however, that because 
the metric is curved, the magnetic field per unit normalized area actually changes and 
becomes weak when |ap gets large. The magnetic field will try to bend trajectories into 
confining circular orbits, but it has to compete with a repulsive potential that tries to 
destabilize the system. Also, the metric has a scaling symmetry where |ap —)■ 7^|ap, 2 : —>■ 7 ^. 
The metric is also positively curved away from |ap = 0, which is a singularity. If we thought 
of the z as an angle coordinate and the |ap as a radial variable, then the curvature wants 
to repel geodesics from hitting |ap = 0. 

The obvious critical point of the equations of motion where nothing moves, is at the 
maximum of the potential and is in unstable equilibrium. This produces a periodic orbit for 




FIG. 8. Parametric plots of {z{9), |ap(0)) for various initial conditions at fixed value of = 0.1. 
We show a BPS trajectory in purple, and examples of scanning over parameters to find periodic 
trajectories in blue and red. The fixed point is marked. The green solution which was found by 
scanning winds twice around the fixed point and it is not BPS. 


any period. Other periodic orbits need to be fonnd by trial and error, and once a snfficiently 
approximate solntion of the periodicity condition is fonnd it is possible to zoom into it. 

Particnlarly simple examples of this search can be performed if the solntions are reflection 
symmetric nnder z —>■ —z. Then the solntion can be characterized by the valne of |ap(0o); 
when z{6q) = 0, and the condition for symmetry forces c^ela^K^o) = 0. We can then move 
the velocity of d 0 z{ 6 o) = ^ as a scanning parameter. Some examples of this procedure are 
depicted in hgure 

When the orbit returns to z{di) = 0 at some later time (not necessarily the hrst time 
around), we can the compute k = de\a\^{ 6 i) from solving the equations of motion, and as 
we scan over ^ we look for sign changes in k. We can then zoom in for the parameter ^ that 
solves the periodicity condition. The point of intersection can be the same one as before or 
it can be different. If it is different we double the time of this half orbit, and it becomes 
periodic. 

It is also interesting to study the hrst order differential equations that arise for states 
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that saturate the BPS inequality at hnite J. These equations are given by 


dez 



dg\a\'^ 



(78) 

(79) 


and are obtained from requiring the vanishing of the sqnares in equation (M after snbstituting 
our ansatz. These can in turn be derived from a variational principle for an auxiliary 
Lagrangian of the form 

L = f-q+l\ogq (80) 

where q = |ap, and z plays the role of the canonical conjugate of q (namely z = 2q = pg). 
The effective one dimensional potential tl{q) = q — ^ log O' is bounded from below and goes 
to inhnity at g —)■ 0 and also at g —)■ oo, so it produces automatically closed periodic orbits 
without self intersections. There is also only one minimum at g = J/4. 

These trajectories are interpreted geometrically as multiply wrapped tori that wrap the 
same torus, where the wrapping number is the number of times we have to go around the 
orbit so that the period matches the number of DO-branes. 

The period of the orbit can be calculated using standard techniques as follows: 


0 = 




(81) 


where the integral is over a periodic orbit, which is characterized by the parametric equation 


w = 9^ + 9 - j log((/) = + q - I log((/) 


(82) 


where W is a constant of integration (the energy associated to the lagrangian function L). 
This solves for z as a function of g readily. Indeed, the parameter W this way dehnes a 
curve with a differential, and then 0 is the period integral over the differential. We have 
to be careful with this interpretation: the curve is real analytic and not a complex curve. 
More general solutions for the BPS states in the continuum limit were found in [26], and 
they similarly show up with various logarithms 

Because U{q) has a non-trivial third derivative about the minimum, it is possible to show 
that the period for orbits very near the hxed point have a decreasing period as we go away 


^ The variable W in their work is related to q in ours, with q ~ 
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FIG. 9. Period of an orbit as a function of |a|max z = 0, and at fixed = 0.1. 

from the fixed point, and for large orbits, the period increases again. This is depicted in 
figure 

Since we need to £x the period to £x the number of DO-branes, while we are allowed to 
vary J, this shows that for some values of J there is more than one non-trivial solution (the 
trivial fixed point orbit can have any period we want to). The larger solution is interpreted 
as a fat torus which in principle can become very large as we decrease j going smoothly to 
the fuzzy spheres in the limit J —)■ 0 while the small one is a thin torus that gets thinner 
and disappears into the trivial solution at hnite j. 

As we increase J, we increase the period of the orbits near the fixed point, as well as 
the minimal period. Eventually the period is too long and the BPS solutions with the hxed 
period we want disappear. This behavior can already be seen for the phase diagram of 3 x 3 
matrices in figure where one of the new family of minima joining transitions i?, C* in figure 
1^ plays the role of the small torus solution whereas the maximal fuzzy sphere family plays 
the role of the large torus. This phenomenon can happen for any multiply wound torus in 
the same way. At the transition C two BPS minima end together by joining with a saddle 
and becoming a non-BPS minimum. Hence, it would not be captured by the BPS solutions 
we are hnding. 


IX. TOPOLOGY CHANGE 

Now that we have the solutions for the matrix ansatz, even if computed numerically, we 
can analyze the geometry of the resulting matrix configurations as membranes using the 
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ideas developed in [19]. The main idea in that work is that when we add a probe brane 
(more precisely, we take a direct sum of the matrix conhguration and a conhguration of 1 x 1 
matrices) and we then consider fermionic off-diagonal modes connecting the probe to the 
matrix conhguration. The resulting Hamiltonian is a truncation of the BFSS matrix model 
to a version with less supersymmetry [2] from which we can obtain a non trivial surface 
embedded in for any set of three Hermitian matrices. The fermions end up carrying a 
non-trivial topological structure that forbids a global splitting into positive and negative 
modes everywhere and thus there have to be degenerations. The geometric locus where 
this is not possible is special in a physical sense; the fermions can not be integrated out. 
When we cross the special locus the number of positive eigenvalues and negative eigenvalues 
of an effective Hamiltonian changes. This spectral how of eigenvalues past zero dehnes 
an integer valued index that colors The interfaces where the index changes dehnes the 
membrane. This is related to the Hanany-Witten ehect m, and more precisely, the creation 
of fundamental strings by crossing D-branes |48j . 

The main calculational tool is to look at the spectrum of the following Hermitian operator 
which depends on the position of the probe characterized by the vector ^ G 

3 




(83) 


2 = 1 


A zero eigenvalue occurs at ^ exactly when det(?^(^)) = 0. At this place, the determinant 
changes sign. We can therefore plot the level set det(H(.^)) = 0 and standard numerical 
algorithms can be used to determine this locus. Since this is a polynomial equation in real 
variables, the corresponding surface is algebraic in nature. In contrast, for the BPS solutions 


at large N, we get a surface that also contains the log function in equation (82). In this 
sense, there is a measurable hnite departure from the finite N and the infinite N limit. We 
can think of this procedure as measuring quantum geometry corrections to large N. The 
most glaring one is that the rotation symmetry group of the solution is reduced from U{1) 
to Zjv, so the matrix solutions are lumpy, and the lumpiness is non-perturbative in iV; we 
get exactly N lumps. 

Our goal stated at the beginning of the paper is to analyze the topology transition from 


a sphere to a torus. As we discussed in section |VHI[ at large N the topology change is 
instantaneous as soon as we turn on a non-zero value of the angular momentum J. As can 
be seen from hgure for the BPS trajectory, this proceeds by forming a very thin funnel 
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between the north pole and the south pole. This funnel represents a condensate of strings, 
as the picture of [TH] suggests we should have. From here we ask a few natural questions. 
The hrst one is if the topology change is apparent as a phase transition in the bosonic set 
of degrees of freedom, that is, if we have to go beyond a phase transition in diagrams 
or[^ The second question is if the transition is instantaneous or not. Lastly we ask what 
happens also when the torus gets very thin (as suggested by the thin ring solution in the 
continuum limit). 


We do not hnd tori for either 2 x 2 or 3 x 3 matrices. This is shown in hgure The 
transition in topology is from a single sphere to three different small spheres. This can be 
understood as a transition from a membrane to a collection of separated DO-branes that 
have been puffed up a little bit from having some off-diagonal excitations, rather than pure 
DO-branes where the ansatz is diagonal. Incidentally, in the trivial conhguration one can 
show that Z commutes with and X~ and they actually commute with each other, so 
this represents a collection of separated DO-branes, because the matrices can be diagonalized 
simultaneously. This is the same type of interpretation as in the BFSS matrix theory |2]. 


The topology transition depicted in figure occurs before any singularity of the phase 
diagram in hgure is encountered, but for the values of J given, it is close to the transition. 


The same procedure for the case of 4 x 4 matrices is depicted in hgure 11 where we see a 



FIG. 10. Transition from a deformed sphere to three small spheres for the case of 3 x 3 matrices. 
The values are at J ~ 1.9 and J ~ 1.99 for the maximal sphere. 


33 










FIG. 11. Transition from a very distorted sphere to a very distorted tours for 4x4 matrices. The 
values are at J ~ 4.24 and J ~ 4.31 for the maximal sphere. 


sphere transitioning to a torns near the phase transition for the maximal sphere, bnt before 

it. 


From these examples it shonld be clear that as far as the bosonic degrees of freedom are 
concerned, the change of topology from a sphere to a torns is smooth. Moreover, it does 
not occnr immediately as in the large N limit. The fermions do detect the topology change. 
The tori that are obtained this way for the maximal sphere for snch small valnes of N are 
very distorted. This shonld improve as we increase N. 

One conld also analyze the geometry and topology in terms of the ideas of fnzzy Riemann 
snrfaces fonnd in [IHlISni- For the analysis of topology, one uses properties of the eigenvalues 
of Z interpreted as a Morse function. This only works for very large N. As far as the 
geometry is concerned, we hnd that the hnite matrices for some sufficiently classical states 


(at sufficiently large N again) would give rise to fuzzy approximations to (82) for some 
values of VF, where q ~ {X~^X~) is interpreted as a matrix and is a normal ordered form of 
the product. These solutions are not algebro-geometric in nature because of the logarithm. 
In this case q and Z commute and can be thought of as classical variables on the torus 
that can be constrained by an equation. Deviations from satisfying the precise equation 
W = + q — l/4Jlog(q') for hxed W should then be interpreted as “quantum corrections” 

in 1/N. These are not quantum in the sense of due to a non-trivial vale of the Planck 
constant h, but should be thought of instead of as quantum corrections to geometry due 
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to the discrete nature of the DO-brane charge. These corrections are due to the small size 
of the M-theory circle in the DLCQ limit [SI]. Here we take the DLCQ of the plane wave 
geometry. 

X. CONCLUSION 

We have analyzed a particularly simple set of periodic classical solutions of the 50(3) 
sector of the BMN matrix model. The solutions are periodic in time modulo gauge trans¬ 
formations. They are also required to preserve the maximal discrete subgroup of rotations 
along the axis of rotation that is allowed by the discreteness of the matrices, namely a Zat. 
At large N, a U{1) of rotations is recovered, so if the solutions have a continuous limit, they 
go to a rotationally invariant conhguration. The solutions are not supersymmetric. As a 
consequence, the rotations that are turned on are not a central charge of the theory. There 
is, however, a BPS inequality that the solutions must satisfy, relating their energy to the 
angular momentum. Some solutions saturate the inequality, and they have a simpler set of 
equations that need to be satished. 

The rigidly rotating solutions can be understood in terms of a system of algebraic equa¬ 
tions. These are in one to one correspondence with critical points of an energy function. 
Solutions found at zero angular momentum correspond not only to the vacua of the ma¬ 
trix model, made of concentric fuzzy spheres, but also other various unstable saddle points. 
These all survive when we turn on the angular momentum. The excitation of angular mo¬ 
tion in the conhgurations always leads to a term in the on-shell energy proportional to the 
angular momentum, rather than starting at the square of angular momentum. This indi¬ 
cates that even though the conhgurations are rotating rigidly, the matrix object can not be 
thought of as a rigid body. 

Both the hnite N and large N results suggest that there is a maximal angular momentum 
beyond which there is only one phase. This should scale like N^. This phase when interpreted 
geometrically corresponds to the brane dissociating into a collection of DO-branes arranged 
symmetrically on a circle. 

The phase diagram of solutions is rather complicated in general and we found it very 
useful to use Morse theory to hnd all the solutions, at least for low values of N. The full 
pattern for a given N has many more saddles than the ones for smaller N. Any solution that 
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is found for lower values of at J = 0 can be combined with other such solutions to build 
solutions at a given N for J = 0, and these are seeds for a family of such solutions at finite 
J. Moreover, for every N there are new solutions, some of which we know, like the maximal 
fuzzy sphere or a maximal fuzzy sphere at half radius. We do not have a systematic way 
to search for the other ones. The parameter space grows in dimension proportional to N 
making it increasingly difficult to find them. Understanding this pattern in general should 
be very interesting. 

At large N, the particular family of solutions we have considered reduce to a variational 
problem for critical points of an energy which is a local integral. The saddle point equations 
reduce to hnding periodic solutions of a pair of coupled second order differential equations, 
while the BPS ones reduce to solving a coupled set of non-linear first order differential 
equations. It is only solutions with the right period that can be used. The period fixes the 
DO-brane charge of the conhguration. 

We also found that the topology change from a sphere to a (sometimes multiply wound) 
torus happens instantaneously in the large N limit, but not so at finite N. Here it is delayed. 
Moreover, the shape of the fuzzy membranes can be very deformed from a circular torus. 
The effects that lead to that deformation are suppressed at large N. It should be interesting 
to investigate this in more detail. The breaking of the symmetry is due to the discretization 
of the DO-brane charge into matrices. In the fuzzy geometry it is a purely classical effect. 
However, in the continuum large N limit, this is supposed to arise from quantum effects 
that are responsible for the quantization of the DO-brane charge. It is often the case that 
quantum effects on D-brane field theories can be captured geometrically, like beta functions 
being captured by brane bending |52l |53] . In this case we find that the discretization that 
appears in matrix theory introduces symmetry breaking effects that are not apparent in the 
continuum limit, and they would not appear in perturbation theory. 

It should also be interesting to understand supersymmetry breaking effects better in these 
solutions and in particular the corrections due to zero point energy. This depends on the 
full theory, as we found that the SO{3) BMN system of classical equations arises in various 
different contexts. Such answers depend on the context. 

Another interesting possibility to examine is that the 50(3) sector of the BMN matrix 
model is also part of the description of the discrete lightcone quantization of the membrane 
in the Penrose limit of AdS^ x 5^ [1], or it’s orbifolds. In particular, one can consider 
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the ABJM model [5l] in the appropriate sector that takes us to the Penrose limit. The 
natural candidates to consider are DO-brane states, which are dual to monopole operators. 
The spectrum of fluctuations around such objects have been analyzed in |S5HSB] and they 
have many fluctuations that are supersymmetric and saturate the BPS bound. It should 
be interesting to try to turn on rotations in AdS^ for such states and see if the extremal 
solutions we have found can be mapped to them at weak coupling as well. This should help 
to understand how the DO-brane theory of the matrix model and the ABJM held theory in 
the presence of monopole states are related to each other. 
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